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Abstract. Recent developments in the physics of low density trapped gases 
make it worthwhile to verify old, well known results that, while plausible, were 
based on perturbation theory and assumptions about pseudopotentials. We 
use and extend recently developed techniques to give a rigorous derivation 
of the asymptotic formula for the ground state energy of a dilute gas of N 
fcrmions interacting with a short-range, positive potential of scattering length 
a. For spin 1/2 fermions, this is E E° + (h? /2m)2TTN pa, where E° is the 
energy of the non-interacting system and p is the density. 



1. Introduction 

Our goal is to find the ground state energy of a low density gas of fermions 
interacting with short range pair potentials. Earlier [^|2lini|4j, the corresponding 
problem for bosons had been solved in the sense that the 'well known' ancient 
formula E/V = 'inp^a in 3D was proved rigorously as an asymptotic formula for 
small p in the thermodynamic limit. Here, p = N/V, where N is the particle 
number, V is the volume and a is the scattering length of the pair potential. 

Owing to the recent interest in low density, cold gases of fermionic atoms, the 
natural question arose: Is a similar formula valid for spin-^ fermions? Of course 
there will be differences, which we can list as follows: 

(1) We can expect that the interaction energy between pairs of fermions of 
the same spin vanishes to leading order in p. The Pauli principle (anti- 
symmetry) implies that the wave function is tiny within the range of the 
pair potential. 

(2) The interaction between pairs of fermions of opposite spin should be the 
same as for bosons. I.e., we expect this energy to be AiraN-^N^/V , where 
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and Ni are the numbers of spin up and spin down particles, respec- 
tively. 

A major difference from the bosonic case is that the energy we are trying 
to compute is only the second order term in the density. The first order 
term is just the kinetic energy of the ideal Fermi gas, which is Eq/V — 

l{6n^f'{iN,/Vr/^ + iN,/Vr/^}. ^ 

Finding an upper bound to the energy in the fermionic case is more dif- 
ficult than in the bosonic case. There, one uses the fact that it is not 
necessary to look for trial states in the bosonic (symmetric) sector be- 
cause the bosonic energy agrees with the absolute minimum energy and, 
consequently, any variational state without symmetry restriction gives an 
upper bound to the bosonic ground state energy. The same is obviously 
not true for fermions, which means that we must rigorously enforce anti- 
symmetry of our trial state. 

The third item gives rise to two problems that did not have to be faced in the 
bosonic case. The first is that we will have to be careful not to 'give up' any low 
momentum kinetic energy in order to control the interaction potential. That kinetic 
energy is part of the main term in the energy and, therefore, has to be left intact. 
Since the interaction range a is very small compared to the Fermi momentum, 
which is proportional to p^^^, there is hope that in the low density regime this 
decomposition of the kinetic energy into low and high density parts can be cleanly 
achieved. This is one of the main technical accomplishments of our work. 

The second main problem is closely related to this kinetic energy consideration. 
It is that the typical particle momentum is of the order of the Fermi momentum, 
which is not zero. Therefore, one could question whether the 'zero energy scattering 
length' contains all the necessary information about the interaction. That it does 
so is a consequence of the good separation of momentum scales mentioned above, 
but proving this fact mathematically requires some thought. 

In the next section we explain the problem and our results in more detail. The 
complete details can be found in and the present note can be regarded as a guide 
to the overall structure of 



2. The Model 

Our model describes fermions interacting through a short range pair potential v. 
The Hamiltonian (in units in which h = 2m = 1) is 



N 








i=l 


l<i<j<N 



The nonnegative pair-potential, ti > 0, is radial (i.e., depends only on \xi — Xj\) and 
has finite range (i.e., has compact support). A favorite physical example is a hard 
core, v{x) = oo for \x\ < a and v = otherwise. 

We generalize the physical problem a little by treating higher spins than i, 
since we let the Hamiltonian Hn act on the fermionic Hilbert space with q spin 
states. Spin i corresponds to q — 2. 
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The particles are restricted to a box A — [0, L]"^ of volume V = L^. The Hilbert 
space describing N fermions with g-spin states in this box is 

N 

n = /\L\k-c''). 

We use Dirichlet boundary conditions. Since the Hamiltonian does not depend on 
the spin degrees of freedom, we may think of the system as consisting of q species 
of spinless fermions. 

We are interested in the ground state energy E{N,L) = inf spec^ in the 
thermodynamic limit, more precisely the energy per unit volume defined by 

eip)= lim 



N/L^=p 

This limit is known to exist by very general arguments dating back to the 50's. 

3. The Main Result 

Theorem 3.1 (Low density asymptotics). For p ^ the minimum energy 
is achieved (to the two leading orders in N/V ) by having equal numbers of particles 
(namely N/q) in each spin state. This energy is 



(3.1) 



where the first term is the energy of a free Fermi gas and the second term depends 
on the pair potential through its scattering length Oy . 

The scattering length of the spherically symmetric potential v is defined 
as the constant such that the spherically symmetric zero-energy scattering solution 

— H — u0 = with lim (j){x) ~ 1 

2 \x\ — *oo 



satisfies 



\x\ 



for all X outside the range of v. 

The proof of Theorem 13. II is given in jS]. 

Remarks: (i) The relevant dimensionless parameter is pa^. Low density means 
that pa^ <C 1 or that p —> 0, when Oy is kept fixed. 

(ii) For a hard core potential Oy is simply the hard core size. 

(iii) From the definition of the scattering length we see, by an integration by 
parts, that 

(3.2) J\\/(t,\' + ^Jv\q^\'=Anay. 

Heuristically, we may interpret this result as saying that for just one pair of particles 
in a large volume the energy is Anav/yol per particle or altogether 87ra^/vol. The 
number of "interacting" fermion pairs is ■^N'^{1 — ^)- The interaction of identical 
fermions can be ignored for low density, as explained earlier. 
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4. The Bosonic Case 



The low density energy asymptotics for bosons (or even without any symmetry 
restrictions, i.e., ii H — ^ L'^{M.^)) was proved by Lieb and Yngvason in [J. 
Some subsequent developments are in [g]-[l4|. For an up-to-date survey of rigorous 
results on the Bose gas j4j is recommended. 

Theorem 4.1 (Dilute Bose gas). For p ^ 



Note that the correlation term is the leading term here. 

In 1957 Dyson |15| proved the asymptotically correct upper bound. This is far 
from trivial since one cannot use a simple product trial function, which would give 
infinity in the hard core case. 

Dyson also gave a lower bound for bosons, which however was not asymptoti- 
cally exact. Dyson's lower bound relied on the following lemma, which, in fact, was 
also a key ingredient in the rigorous lower bound in pp. 

Lemma 4.2 (Dyson's Lemma). Assume v supported in {\x\ < Rq}. Let 9ji 
be the characteristic function of {\x\ < R}. Then for any positive radial function 
U supported in Rq < \x\ < R with J U = At: we have the operator inequality 



Dyson's Lemma states that one may, as a lower bound, replace w by a much 
softer potential w = OyU with the property that (47r)~^ J — (4''')^^ J o-vU — a^. 
However, one has to use all of the kinetic energy to do this. 

In the rigorous argument in 1 only almost all the kinetic energy is used with 
the Dyson Lemma. The rest is used to carry through a perturbation argument. 

As explained in the introduction we cannot afford to give up even a fraction 
of the low momentum part of the kinetic energy in the fermionic case, since it is 
responsible for the leading free Fermi term. We therefore need a generalization of 
Dyson's lemma, which we describe below. 

Early History of the Bosonic Problem. Lenz |16| was the first to calculate 
the energy H4.1|l . using the heuristic pair argument given in Remark (iii) after 
Theorem 13.11 

Bogoliubov |17| did a perturbative expansion, but got (Stt)^^ J v instead of 
fli,. He realized, of course, that this was meaningless (e.g. for the hard core). As 
Bogoliubov points out in his paper. Landau realized the connection to the scattering 
length. Namely, that {8tt)~^ J V is the first term in the Born series for a„. See also 



Using Bogoliubov's method, the next term in the expansion beyond Airavp^ has 
been calculated to be 



An important open problem is to prove this formula rigorously. As in the leading 
term, this correction depends on v only through the scattering length. 



(4.1) 



e^{p) =Ana^p^ + o{p^), 




m 
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5. The Generalized Dyson Lemma 

The sohition to the problem that we cannot give up any fraction of the low 
momentum kinetic energy is to give up almost all the energy in the large momentum 
regime only. Large here means momenta much greater than the Fermi momentum. 
If we give up kinetic energy corresponding to momenta greater than the Fermi 
momentum only, and in such a way that as a function of momentum the kinetic 
energy is still monotone, we do not affect the leading free Fermi gas term. The 
following lemma, which allows us to do this, is the main technical advance in our 
work. 

Lemma 5.1 (Generalized Dyson Lemma). Assume v to be supported in 
{\x\ < Rq}. Let 6ji be the characteristic function of {\x\ < R}. Let x be a radial 
function, < x{p) ^ 1^ such that the Fourier transform h{x) = 1 — x{^) bounded 
and integrable. Define 

= sup \h{x ~ y) ~ h{x)\ and wb{x) ^ \fR{x) [ fR{y)dy. 
\y\<R Jr3 

Then for any e > and any positive radial function U supported in Rq < \x\ < R 
with J U = Att we have the operator inequality 

-^x{p)Or{x)x{p)^ + H^) > (1 - e)a,U{x) - ^wr{x). 

i? is a parameter we shall choose such that <$:. R <$:. p^^^^. We will choose the 
function x to be approximately the characteristic function for the set of momenta 
greater than the Fermi momentum, which is proportional to p^^^. This corresponds 
to decomposing the kinetic energy into low and high momentum parts as shown in 
the figures below. 



Low momentum regime 



p'^xip) 

High momentum regime 




6. A Lower Bound to the Energy 

To explain how the generalized Dyson Lemma is used to get a lower bound to 
the energy we consider, for simplicity, the case q = 2 and write X = {xi, X2 ■ ■ ■ ,xn-^), 
Y — (yij 2/2 ■ • • VN^) for the coordinates of the two species of fermions. Without loss 
of generality we may assume that = iV|. The reason is that the system is 
SU{2) symmetric and every state of total spin S can be rotated to a state of total 



6 



ELLIOTT H. LIEB, ROBERT SEIRINGER, AND JAN PHILIP SOLOVEJ 



Sz ~ 0. Nevertheless we will keep the particle numbers arbitrary for the moment 
for pedagogical clarity. 

We first consider the Y particles as fixed. For a lower bound we may ignore the 
repulsion between like fermions since v > 0. From the generalized Dyson Lemma 
we obtain the following result. 

Corollary 6.1 (Generalized Dyson Lemma in multi-centered case). 

Let R > be bigger than the range of the potential Rq. If \yi ^ yj\ > 2R for all 
i j, then 

-Vx(p)^V + 5 ^ w(a; - > ^ (^{1 - e)ayU{x - yi) - —wr{x - yi)^ . 

i=l i=l ^ 

To use this we must ignore y-particles that are closer than 2R apart. This can 
be controlled by the following a-priori bound. 

Lemma 6.2. Let Ir{Y) be the number ofy's that are less than 2R > from 
their nearest neighbor. Then for any fermionic wave function "^{Y) 

/i?^') < C•^(^')i?^ where r(*) = ^(V^,*, V^,*) 

i 

and C is a constant. 

This lemma is an easy consequence of a theorem jl9l Theorem 5] relating 
nearest neighbor distance to kinetic energy. 

For an approximate ground state \E' we can assume 
that we wiU choose Rp^/^ < 1, and thus {'1',Ir(Y)^) < TV. 

The soft potential can now be studied with the remaining low momentum 
kinetic energy. 

[-^^(1 - xb))V., + ((1 - e)ayU{x, - y,) - ^WRix, - t/,))] , 

2 — 1 

where refers to the sum over the Ni — IniY) y-particles that are a distance 
greater than 2R from their nearest neighbor. Since we can ignore Ir{Y) the sum 
of the lowest eigenvalues of the operator in [ ] is approximately 

^ (6^^)'/' iVt^/^L-^ + N^^N^-L-'a, J U . 

We get a similar contribution when we consider the Y particles as movable and the 
X particles as fixed. If we also use that J C/ = 47r the sum of the two contributions 
will be 

Here we see again that the minimum occurs for = , and is given by (|3.1|) . 

7. The Upper Bound 

For an upper bound we first construct a Bijl-Dingle-Jastrow function for par- 
ticles of the same spin: 

G{X)^\{g{x.,-x,), 
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where g{x) = if |a;| < Rq, g{x) = 1 if |a;| > s, and i?o ^ s *C p^^/^. The role of 
G is to ensure that hke fermions do not interact. 

Next, we construct a Bijl-Dingle- Jastrow function for particles of opposite spin: 

(7.1) F{X,Y)^l[f{x.-y,), 

where / (which is not normalized, since it, like g, satisfies f{x) — 1 as |a;| — > oo) is 
chosen so that 

(Compare with (|3.2(l .) The role of F is to give the correct pair energy between 
different fermions. 

Finally, the trial state we shall use is 

*(X, Y) D{X)D{Y)G{X)G{Y)F{X, Y) , 

where D = u^j^ A • • • A is a Slater determinant of eigenfunctions Ua of the 
Dirichlet Laplacian. The choice of this 4" is quite obvious. The problem is that it 
is not normalized. The normalization problem can be addressed with the aid of the 
following combinatorial lemma, for which we claim no originality. 

Lemma 7.1 (Key combinatorial lemma). Let $ = 0i A • • • A denote a 
Slater determinant of n linearly independent functions (j)aix). Let M denote the 
n X n matrix 



(7.2) Map - J (j)l{x)(j)p{x)dx. 

(i) The norm of ^ is given by = detM. 

(ii) For 1 < m < n, the normalized m-particle density of $ is given by 

^ ^ [0(a;)] A • • • A [(p{x^)] 



{xi,...,Xm)^ <^[(l)ixi)]A--- A[(j){Xm)] 

Here, [cf){xi)] denotes the n-dimensional vector with components 4>a{x) ■ 

We use this lemma, for fixed Y , with (t>a{x) = Uaix) Y[j f{x — yj). The matrix 
M in H7.2|l then depends on Y and is denoted by M.{Y). The following estimate 
shows that although the 4)a are not orthonormal, M.{Y) is close to / if the separation 
between the F-particles is not too small. 

Lemma 7.2 (Key estimate). Assume that \yi — yj\ > s for all i ^ j. Then 
\\L — M{Y)\\ as s/oy oo and N^^'^s/L 0, uniformly in N and Y . 

Note that if / were identically 1, then A1 would be equal to the identity. 
The proof of Lemma 17.21 uses the Poincare inequality to estimate the effect of the 
deviation of / from 1 on the integral in (|7.2|l . 

The assumption that \yi — yj\ > s in Lemma 17.21 is satisfied since our trial 
function vanishes otherwise, thanks to the property of G. Lemma 17.21 is concerned 
only with the effect of F on the norm of 4". The effect of G on the norm of 
has to be controlled by different means. For this purpose we find it necessary to 
break space up into boxes (whose size is large but independent of L) and confine 
the particles to these boxes, i.e., regard particles in different boxes as independent. 
We refer to [5] for details. 
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